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Abstract 

A simple introduction to physics of CMBR polarization and the Faraday rotation 
of the latter in cosmic magnetic field is presented. The content of the lecture is the 
following: 1. Description of polarization of photons. 2. Polarization field of CMBR. 
3. Faraday effect. 4. Cosmic magnetic fields. 5. Faraday rotation of CMBR 
polarization. 
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1 Description of polarization of photons 

The polarization state of an electromagnetic wave is determined by the vector of its electric 
held which is orthogonal to the direction of the propagation. The polarization density 
matrix is dehned as the average value of bilinear combinations of the components of the 
electric held of the propagating wave [T]: 

Pii = (£,£:•> (1) 

Evidently, pij is the 2nd rank tensor in 2-dimensional {x, ?/)-space if photon propagates 
along z. This tensor has the following well known invariants: 

1. Trace, which is equal to intensity of radiation: 

T = SijPij = \Ex\'^ + ( 2 ) 

2. Helicity: 

V = €ijPij (3) 

Non-zero helicity V corresponds to circularly polarized photons when the average value of 

the photon spin projection on the direction of the photon propagation is non-vanishing. 

As any 2 x 2-matrix pij can be decomposed in terms of the unit matrix, I and three 
Pauli matrices, cr^, k = 1,2, 3. The coefficients in front of the latter are called the 
Stokes parameters: 


Pij — T (//2 -|- ^fcCTfc) (4) 

For a particular and practically interesting case of photons with a hxed frequency the 
components of electric fields can be written as: 

= Eqc^ exp[-iujt + if3^] 

Ey = EoCy exp[-iut -f i/3y], (5) 

where the components of the photon polarization vector satisfy the normalization 
condition -|- = 1. 

The Stokes parameters for this case can be expressed as: 

6 = e^ey sin(^,j, - py) (6) 

where ^2 is invariant with respect to rotations and pseudoscalar with respect to mirror 
reflection. It describes circular polarization, i.e. the photon helicity: 


A = sk/o; 

The other two Stokes parameters .^1 and .^3 describe linear polarization: 

6 = (ex-4)/2 

= e^Cy COs{Pa: - Py) ( 7 ) 
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They transform under coordinate rotation in (x, ?/)-plane by angle 0 as: 

= .^1 cos 20 — ^3 sin 20 

^3 = 6 sin 20^3 cos 20^3 (8) 

Making proper rotation one can always arrange vanishing of one of these Stokes parame¬ 
ters, e.g. 

= 0 ( 9 ) 

Using the transformation law (jH)) one can check that the following combinations of the 
Stokes parameters are the eigen-vectors of rotation: 

^3 ± iii exp [±2i0] (^3 ±i^i) (10) 


2 Polarization by Thomson scattering 


As is well known, see e.g. the book P], elastic scattering of unpolarized photons on 
unpolarized electrons 

'y + e ^'y'+ e' (11) 

produces polarized photons. The reaction amplitude can be written as 

A = e'Ai, (12) 

where e' is the polarization vector of photons in the final state. The polarization matrix 
of the scattered photons, up to normalization factor, is expressed as: 


Pij -i-' -j 


AiA* ( 13 ) 

If we choose the coordinate z in the direction of 7' and x in the reaction plane and denote 
by 9 the scattering angle then the only non-zero Stokes parameter would be equal to 

sin^ 9 sin^ 9 


^3 = 


( 14 ) 


a;/a;'-I-a;Ya; — sin^ 0 1-|-cos^ 0 

where uj and oj' are the energies of the initial and hnal photons respectively and the 
approximation of non-relativistic electrons is made. 

The result m can be easily understood. The only non-vanishing combination involv¬ 
ing e' in the amplitude is: 

e' k ~ sin 9 ( 15 ) 

Hence, from eq. m follows that 

^3 ~ (sin 6*)^ ( 16 ) 

The other Stokes parameter, which may be non-zero, vanishes by the choice of coordinate 
direction, ^1 = 0, while ^2 = 0 due to parity conservation. 

The degree of polarization of the scattered radiation is proportional to the differential 
Thomson cross-section, which is equal to P: 

da 3 87ra^ 
dO, IGtt 3 m^ 


(1 -I- cos^ 9 ) 


(17) 
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3 Polarization of CMBR 


According to the previous section the photons of CMBR could become polarized due to 
Thomson scattering on the electrons in cosmic plasma. However, the polarization must 
vanish in homogeneous and isotropic world. It is evidently true, simply because there is 
no preferred direction in such a world. In realistic cosmology small but non-zero density 
perturbations existed in the primeval plasma prior to hydrogen recombination. As we 
see below, the polarization of CMBR would be induced by the the inhomogeneities in 
the distribution of photons, i.e. by the angular fluctuations of their temperature. This 
effect was noticed many years ago, the pioneering papers include refs. |2]. For review and 
history see e.g. [3]. 

We assume that circular polarization vanishes. As we mentioned above, see also ref. j3], 
such polarization cannot be induced due to parity conservation. Bounds on possible 
new interactions violating parity from the absence of circular polarization of CMB are 
discussed in ref. |S] . With the present day accuracy they are very far from the usual weak 
interactions. 

As we have seen above, in the usual case of the linearly polarized CMBR radiation 
held the intensity of polarization is described by two functions, Q = and U = T^i, 
where T is the total intensity of radiation with frequency u: 


P 


Q u 

-U Q 


The total polarization originated as a result of the Thomson scattering should be 
obtained by the integration of expression (HI with account of (HU over all angles dQ = 
d cos 9d(j) with rotation around z to the common coordinate system: 

Q — iU =— [ do; sin^ 0 exp[2i0] T (6*, 0) (18) 

o'n J 

where ctat is some normalization area over which the data are averaged. Thus one can see 
that the polarization of the scattered radiation is proportional to quadrupole moment of 
the incoming radiation. 


4 Properties of CMBR polarization field 

For description of the CMBR polarization is convenient to know that in addition to 
the invariants o and (0 there exist two more invariants which can be constructed by 
differentiating of the polarization matrix. One of them is a scalar: 


S = didjpij, (19) 

while the other is a pseudoscalar: 

P = CikdidjPjk (20) 

If density perturbations are created by a purely scalar field, then the only way to 
construct a vector out of this scalar is to apply the operator of differentiation, i.e. gradient. 
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dj. A second rank tensor can be either obtained by the second order derivatives, dj dt or 
by multiplication by the invariant Kronecker tensor 6ij. In particular, the traceless part 
of the polarization tensor should have the form: 

pij = {2didj - ( 21 ) 

where \1/ is a scalar function of coordinates. 

One can easily check that for the case of scalar perturbations the pseudoscalar invariant 
(P|) vanishes: 

P = 0 (22) 

Thus we arrive to an important conclusion that a non-zero P is an indication for tensor 
perturbations or, more generally, for something beyond scalar perturbations. A short list 
of those includes: 

1. Vector perturbations, e.g. magnetic helds. The relevant part of the polarization 
matrix can be written in the form: 

pij = diVj + djVi (23) 

where V describes a vector perturbation held. Correspondingly, P may be non¬ 
vanishing and equal to: 

Pv = eijd'^diVj (24) 

2. Tensor perturbations, which might be produced by gravitational waves. In this case 
the CMBR polarization matrix may contain the contribution: 

Pij dih^j -|- djh^ii ( 2 ^) 

and the pseudoscalar invariant takes the form: 

Pt = eijd'^dihsj (26) 

3. The pseudoscalar invariant may be non-zero even in the case of purely scalar per¬ 
turbations but in the second or higher orders pin. The polarization matrix would 
be proportional to the products of gradients of two different scalar functions and 
^ 2 : 

p,, = d,^,dj^2 + di'^2d,<il,, (27) 

where e.g. T 2 = The P-invariant is equal to: 

P 2 = eij {did ‘^^1 a,T 2 + did ‘^'^2 5 ,^ 1 ) (28) 

The decomposition of the polarization matrix according to the magnitude of the scalars 
S and P is considered in ref. jS]. It is analogous to magnetic/electric or B/E - decomposi¬ 
tion P but the former looks more natural because no vector held of electric or magnetic 
type is associated with 2x2 polarization matrix. 

However, there may be other, non-primordial, sources for non-zero P due to propaga¬ 
tion of CMB photons through the large scale structure. They are analyzed in ref. nm. 
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5 Eigenvectors of polarization matrix 

In some works polarization fields is characterized by the “vector” 

W = {Q,U) 

However, W is not a vector bnt some mixtnre of components of a second rank tensor. 
Thongh snch a description is possible, it is not convenient HU. Field lines of such “vector” 
change under rotation in a non-trivial way. We see portraits of famous astronomers and 
the portrait of Louis XIV in this beautiful room. If someone takes one or other portrait 
and rotates it, or ever simpler, just turns one’s own head, the picture would surely remain 
the same. This would not be true for the pictures of CMBR polarization held presented 
in terms of “vector” W. 

Much more convenient is to present polarization maps in terms of held lines of the 
eigenvectors N of the polarization matrix pij H21I5- interesting feature of such maps 
is unusual singularity points of the held lines of the eigenvectors N. As is known from 
the classical analytical mechanics, dynamical systems may normally have the following 
three singularity types: saddle, focus, and knot. However, the eigenvectors N are not 
analytic near singularity, where |N| = 0. Because of that the character of singularities 
changes. There still remain three types of them but now they look diherently, see hg. 1. 
We suggested the following names for them: saddle (the same as above but now it has 
three, instead of four, straight line asymptotic), beak, and comet H21E1- analysis of 
statistical properties of these singularities is performed in the papers El HU 112 ESI- 

In hg. 2 the simulated map of CMBR polarization is presented, according to ref. El- 
One can easily see the singular points corresponding to the three types indicated in hg. 
1. In contrast to description of the direction of the polarization held in terms of W, the 
map of held lines of N remains the same independently of the angle at which one observes 
it. On the other hand, as we mentioned above, diherent singular points of the map of the 
held lines of W transform into each other under rotation and the general picture becomes 
completely diherent. 

6 Faraday effect 

The Faraday ehect is a rotation of the polarization plane of linearly polarized photons 
in magnetic held B, due to interaction with electrons in the medium. Such rotation is 
generated by the breaking of the rehection symmetry in presence of magnetic held. In a 
medium without rehection symmetry refraction indices for left- and right-handed photons 
are diherent, n+ ^ ri-. Linearly polarized electromagnetic wave can be decomposed 
into two rotationally polarized ones rotating in opposite directions. We assume that the 
electric held of the initial wave was directed along axis x. Then using the simple identity: 

1 = (l + i)/2 + (l-i)/2 (29) 

we expand in terms of two helicity eigenstates of the incoming photon: 

= E^ = {E+ + E_)/2 (30) 
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a b c 

Figure 1: Flux lines for three different types of singular points: (a) saddle, (b) beak, and 
(c) comet. Dashed lines show peculiar solutions, separatrix. 



Figure 2: Simulated map of CMB polarization eigen-vector held N Solid lines show the 
hux line behavior near singular points where polarization vanishes. 
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Each helicity state propagates independently and acquires phase proportional to the dis¬ 
tance of propagation: 

= exp {ik±l) (31) 

Since the refraction index for different helicity states is different by assumption, then 
k± = ko ± Ak and 

exp[ikol] cos{Akl) 

Ef^') = exp[iA:o/] sin(AA:/) 

The relative phase remains zero. It means that the initially linearly polarized wave remains 
linearly polarized but the polarization plane is rotated by the rotation angle: 

$ = arctan = ^^1 

Thus the calculation of rotation angle is reduced to the calculation of the refraction 
index of ionized gas - this is the cosmic medium where photons of CMBR propagate. 
The equation of motion of electrons, with charge (—e), in superimposed external (large) 
magnetic held Bq and (weak) electromagnetic wave Eexp(ici;t) reads: 

f = eBo X f — eEexp[ici;t] (32) 

To solve this equation we need to decompose the propagating wave in terms of helicity 
states: 


E = ^7+ (n. 




C'_(n. 


iw„ 


for which the equation diagonalizes and can be solved as 

eE-t 


x± = 


muj {u =F <x’b) 


where ojb = eBo/m. 

The electric polarization moment of the plasma is easily found: 


(33) 


(34) 


V± = -Neex± 


(35) 


and correspondingly the dielectric constant is 


e-t = 1 -I- 4:7iV/E 


1 + 


dvre^iVe 
mu {u =F ^b) 


The plasma refraction index for different helicity states of photons is given by the 
standard expression: n± = y/e± and thus the differential Faraday rotation angle is equal 
to: 


d(j) 27rNee^Bo 

dl 


where m is the electron mass, e'^ = a = 1/137, and Ng is the number density of electrons. 
Usually the results is presented in terms of frequency u = u/{2'k) or wave length A = 1/v. 





7 Cosmic magnetic fields 

The origin of large scale cosmic magnetic fields remains one of cosmological mysteries. 
They surely exist in galaxies with the coherence scale of a few kpc and rather large held 
strength: 

Bgai = a few /iG, (37) 

for a review see ^3], more recent data can be found e.g. in ref. ng. Less certain are the 
observational data in favor of existence of intergalactic magnetic helds at the scales ~ 1 
Mpc. Still there are rather convincing indications in favor of the latter with the amplitude 
Big ~ lQ~^Bgai. It is rather interesting that the amplitude of these helds is related to the 
galactic held by the inverse ratio of the corresponding scales squared. This would be so if 
the adiabatic compression took place and the held was amplihed as B l/t^. Since the 
ratios of scales are 



one would expect Bgai ~ if galactic and intergalactic magnetic helds have a 

common origin and galactic dynamo amplihcation was not ehective. The latter might be 
quite efhcient. According to diherent estimates CB] the amplihcation factor could be as 
large as If galactic dynamo indeed amplihed primordial seed helds by this large 

factor, it means that the original helds were too weak to be observable by the Faraday 
rotation of the CMBR polarization. Otherwise, if Big ~ 10“® Gauss, the ehect may be 
noticeable. 

There is rather strong evidence accumulated recently im in favor of quite strong 
magnetic helds in galactic clusters with the magnitude close to those in galaxies. If this is 
indeed true, then the chances of observation of the Faraday rotation of GMBR polarization 
would be even higher. 

A vast literature exists on possible mechanisms of generation of primordial magnetic 
helds (for a review see e.g. ^H])- Diherent mechanisms can be roughly speaking separated 
into the following classes: 

1. Galactic processes, stellar phenomena and reconnection of held lines. 

2. Processes during structure formation. 

3. Processes during the recombination epoch; vorticity, VxD, may be generated in 
the second order in density perturbations. 

4. Processes in the early universe: 

a) inhation, might produce small helds but at large scales; 

b) phase transitions could create large helds but at small scales. 

All these mechanisms, except for the hrst one, might create noticeable helds at the 
GMBR decoupling which can be potentially observable by the Faraday rotation of the 
GMBR polarization. This ehect was hrst discussed in ref. and attracted considerable 
attention 120112111221 during the last years. 
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8 An estimate of the rotation angle 

The dominant contribution to the rotation angle comes from the period near the hydrogen 
recombination epoch. Indeed, the dependence of the rotation angle on the cosmic scale 
factor is: 


~ }?NeB adT] a? 


1 1 


( 39 ) 


where rj is the conformal time, the magnetic held is assumed to evolve as B ~ 1/a^ 
(adiabatic compression and no dynamo amplihcation). Thus the rotation is dominated 
by an early epoch. Before the recombination the photon mean free path along which the 
rotation angle could be “accumulated”, I free, is very small and practically (<h) = 0. After, 
the recombination the number density of free electrons, N^, drops down by the factor 
10“^. Thus the most favorable period for the rotation of the polarization plane in large 
scale magnetic helds is the recombination epoch. 

The differential rotation angle is equal to: 


dr] 


XeNeC^a 

271711 ? 


Bn 


( 40 ) 


where x^. is ionization fraction and n is the unit vector in the direction of propagation of 
radiation. By assumption of adiabatic compression Ba^ = const = BqoI, i.e. it is equal 
to the present day value. 

The optical depth is expressed through the Thomson cross-section ut and the number 
density of free electrons and is equal to 


dr 

dr] 


= NedTO 


(41) 


The total rotation angle (for homogeneous held along photon propagation) is easily esti¬ 
mated: 


$ = 


SAqBo ■ n /■ , , , 3AnBn ■ n 


Ibvr^e 


dr exp (—r) = 


'0-^0 


Ibvr^e 


( 42 ) 


(here the usual in this held, but unusual for particle physics, convention = a is used). 
Numerically we hnd: 


10“®Gauss J \ no ) 


(43) 


(it helps to know that 1 Gauss = 6.9 ■ 10 MeV^) 


9 Statistical properties of magnetic fields and rota¬ 
tion angle 

We express all the quantities below in terms of the present day values, in particular: 

Bo (x) = 0 ^( 77 ) B(x,r 7 ). (44) 
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As usually the field amplitude is expanded in terms of its Fourier modes: 

Bo(x) = -^ / 

Below we omit sub-zero for simplicity of notations. 

The correlator of the field amplitudes is: 


{Bi(xi)B^(x2)> = Ci^(|xi - X 2 I) 


( 46 ) 


The correlator depends on the difference of the coordinates because of the average ho¬ 
mogeneity and isotropy. Correspondingly the Fourier modes should be delta-correlated. 
Their correlator can be expressed through two scalar functions A{k) and S{k) |2dj : 

n,, = (6,(ki)fo-(ka)) = 2 {27rf6 (ki - k2) 

[(5p - KiKj) S{k) -F ieijiKiA{k)] (46) 


where Ki = fci/|k|. 

The energy of magnetic field depends only on S{k): 

j d^xB^ ~ ~2 j dkk‘^S{k) (47) 

while A{k) determines the so called helicity of magnetic field 

The correlator of the rotation angle can be expressed through the correlator of mag¬ 
netic field as: 

($(n)$(m)) = 

j dr]'g{r]'){[Bo{Ar] n) ■ n)(Bo(Ar7'm) ■ m]) 

where g{r]) = {dr/dr]) exp[— r(r7)] and dg = g — g^. 

Using eq. (HU|) one finds 

((B • n)(B ■ m)) = j d^k {[{nm) - (n ■ K)(m ■ k)] S{k) + 

i [(n X m) ■ k] A{k)} exp [—-ik (nAr; — mAr;')] (48) 

It can be checked that the term containing A vanishes (SUES- It prevents from measure¬ 
ments of helicity of magnetic field through Faraday rotation, in contrast to the results of 
earlier papers 

Performing the necessary integration we obtain: 


((Bn)(Bm)) = [(nm)C^(r) -|- (nr/r)(mr/r)(C||(r) — Cj_(r))] , 
where r = nAr; — mAg', and 


(49) 


CAr) = 


3(27r)3 


dk£^{k) 


jo{kr) 


j2{kr) 


CAr) = 


3(27r) 


dkE^ik) [io{kr) + j 2 {kr )], 


(50) 
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Here ji(x) are the spherical Bessel functions of the order and S^ik) is the magnetic 
power spectrum; 


^B = J nftyodk£B{k), 


SB{k) = 27rk^S{k) 

(51) 

The function S{k) is usually parametrized as 


S{k) = Sok^^ exp[-{k/K)^] 

(52) 


For this simple form of S{k) the integral can be taken analytically. The correlation 
function of the Faraday rotation measure, according to ref. [21] is presented in fig. 3. The 
results are not particularly sensitive to the power ns- 

RR’(0) 

RR’(O) 

1 

0.8 
0.6 
0.4 
0.2 
0 

- 0.2 


Figure 3: Faraday rotation measure correlation RR'{6) as a function of the separation 
angle 6. The three lines correspond to the magnetic field spectral index ns = 2 {solid line), 
ns = 4: {dashed line) and ns = Q {dotted line). The correlation length of the magnetic 
held is taken as ^ = 20 Mpc. 
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10 Conclusion 

1. We do not understand how large scale cosmic magnetic helds have been formed. If 
Bgai and Big have the same origin and galactic dynamo did not operate, impact of 
primordial helds could be observable in CMBR polarization. Another way around. 
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an observation of the Faraday rotation of the CMBR polarization can bring an im¬ 
portant information about primordial magnetic helds and help to solve the problem 
of their origin. 

2. P-type (or P-type) polarization created by magnetic helds may mimic gravitational 
waves but they could be possibly distinguished due to different frequency depen¬ 
dence. 

3. Eigenvector description of the polarization held may be useful. Their statistical 
properties may depend upon the strength of the primordial magnetic helds. This 
problem deserves further consideration. 
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